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Abstract: We examine the behavior of the retarded Green’s function in theories with Lif¬ 
shitz scaling symmetry, both through dual gravitational models and a direct field theory 
approach. In contrast with the case of a relativistic CFT, where the Green’s function is fixed 
(up to normalization) by symmetry, the generic Lifshitz Green’s function can a priori depend 
on an arbitrary function G{uj), where ui = uj/\k\^ is the scale-invariant ratio of frequency to 
wavenumber, with dynamical exponent c. Nevertheless, we demonstrate that the imaginary 
part of the retarded Green’s function (i.e. the spectral function) of scalar operators is ex¬ 
ponentially suppressed in a window of frequencies near zero. This behavior is universal in 
all Lifshitz theories without additional constraining symmetries. On the gravity side, this 
result is robust against higher derivative corrections, while on the field theory side we present 
two z = 2 examples where the exponential suppression arises from summing the perturbative 
expansion to infinite order. 



Contents 


1 Introduction 1 

2 The Green’s function in a scale invariant theory 4 

3 Holographic Lifshitz models 5 

3.1 Bulk higher derivatives 7 

3.2 Consistency of the higher derivative expansion 8 

4 WKB analysis of the spectral fnnction 10 

5 Field theory models with z = 2 17 

5.1 The quadratic band crossing model 18 

5.1.1 Renormalization group analysis 19 

5.1.2 Boson correlation functions 19 

5.1.3 Dirac theory revisited and systems with higher 2 ; 22 

5.1.4 Limitations of the analysis 23 

5.2 The quantum Lifshitz model 24 

6 Discussion 26 

A Perturbative expansion of the WKB integral 28 

B The free quadratic band crossing theory 30 


1 Introduction 

The AdS/CFT correspondence provides us with a remarkable strong-weak coupling duality 
between a bulk gravitational theory and a boundary field theory in one fewer dimension, and 
as such it has found numerous applications in calculating interesting observables of strongly 
coupled field theories using weakly coupled holographic methods. In recent years, the number 
of proposed gravity duals to interesting strongly coupled QFTs has increased significantly. In 
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addition to the relativistic case of asymptotically AdS backgrounds, there exist interesting 
gravitational duals exhibiting non-relativistic (or Lifshitz-) scaling symmetry [1], as well as 
Schrodinger symmetry [2-4]. These non-relativistic backgrounds are particularly relevant 
for studying condensed matter systems at strong coupling, where exact analytic results are 
difficult to obtain using traditional methods. 

In many standard cases, the mapping between bulk and boundary is easily obtained. 
This is especially true when there is a brane interpretation, such as the familiar picture of 
IIB theory on AdSs x S^. In this case, the AdSs supergroup SU(2,2|4) is identical to the 
superconformal symmetry group of the four-dimensional AA = 4 super-Yang Mills theory. As 
a result, all observables are constrained by the superconformal symmetry, and in particular 
the two-point functions are fully determined up to normalization. For example, the retarded 
scalar Green’s function in momentum space must have the form 

GR{q^)=A{-q^)^-\ q^=u^-\k\\ ( 1 . 1 ) 

where A is an overall constant and A is the conformal dimension of the scalar operator 0/\. 

The mapping between condensed matter systems and backgrounds with non-relativistic 
scaling symmetry is often less obvious. In this case, we must often fall back to the general 
strategy of constructing a holographic dual to a given field theory by matching symmetries and 
conserved quantities [1-4]. Moreover, non-relativistic scale invariance is no longer sufficient 
to fully constrain the form of the two-point functions. Consider, for example, the case of 
Lifshitz scaling with dynamical exponent z, where energy and momentum scale as a; —)• 
and k —)> Afc, respectively. This scaling symmetry only constrains the form of the Green’s 
function up to an arbitrary function of the scale-invariant quantity Cj = a;/|A: |^: 

Gjj(u;,fc) = |fc|2-g((h). (1.2) 

Here v is the energy scaling dimension, and the momentum-dependent prefactor is chosen to 
give Gr the proper scaling dimension. 

The form of the Green’s function (1.2) holds for any (isotropic) scale-invariant theory, 
whether computed directly from the field theory or via the holographic dual. However, in 
general, Q{Cj) cannot be fixed by matching symmetries alone. (If additional symmetries are 
imposed, such as z = 2 Schrodinger symmetry, then the Green’s function may become fully 
determined.) This suggests that symmetries are not sufficient for connecting non-relativistic 
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theories to their holographic duals, and in particular that the duality map must include 
additional dynamical information. 

At the same time, the bulk theory yields a preferred choice of the Green’s function 
obtained from the classical two-derivative bulk action. For z = 2 Lifshitz, the holographic 
scalar Green’s function was obtained analytically in [1] , while a WKB calculation for arbitrary 
z > 1 demonstrated a characteristic exponential suppression of the spectral weight (i.e. the 
imaginary part of the Green’s function) in the limit ci) —)• 0 [5]. This has been interpreted 
as an “insensitivity” of the boundary theory to small changes of the geometry near the 
horizon. The same exponential behavior is responsible for making the smearing function of 
both Schwarzschild-AdS and Lifshitz spacetime a distribution rather than a true function 
[6, 7], and has been interpreted as a loss of bulk locality for such non-relativistic geometries 
[ 8 ]. 

It is natural to expect that different field theoretic models with the same dynamical ex¬ 
ponent z will yield different Green’s functions. This raises the issue as to how the holographic 
dual can distinguish among these models. For unbroken scaling symmetry, the bulk geometry 
is essentially fixed to be pure Lifshitz (if we work within the context of general relativity; see 
e.g. [9-11] for other approaches). Thus the background alone cannot distinguish between dif¬ 
ferent models, and we are mainly left with the dynamics of the bulk fields as the distinguishing 
characteristic. In particular, the addition of higher derivative terms to the bulk equations of 
motion will directly affect the form of the holographic Green’s function. This is in contrast 
with the relativistic case, where higher derivative corrections may affect the constant A in 
(1.1), but will not otherwise modify the functional form of the retarded Green’s function. 

Once we allow for a higher derivative expansion in the bulk, it may seem that some 
predictive power is lost, since the holographic Green’s function would in principle be sensitive 
to all of the infinitely many higher derivative terms. However, we demonstrate that there are 
universal features that remain. In particular, the characteristic exponential suppression of 
the spectral function in the low frequency regime found in [5] is robust with respect to higher 
derivatives in the bulk, as long as the frequencies stay above a (momentum-dependent) cutoff. 

We furthermore show that this exponential suppression arises in field theory models with 
z = 2 scaling. In particular, for both the quadratic band crossing model of [12] and the 
quantum Lifshitz model [13], a simple kinematical argument demonstrates that the exponen¬ 
tial suppression arises because one has to go to higher and higher orders in the perturbative 
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expansion to see non-zero spectral weight in the limit w —)> 0. 

This paper is organized as follows. In section 2, we briefly review the form of the scalar 
Green’s function in a theory with Lifshitz scaling. Then in section 3, we set up the compu¬ 
tation of the Green’s function in a holographic Lifshitz model with bulk higher derivatives. 
Following that, we perform a WKB analysis of the spectral function in section 4 and show 
that it has a universal exponential suppression at small a). In section 5, we demonstrate that 
this exponential suppression can be seen directly from a held theory perspective. We focus 
on the quadratic band crossing model, but also consider the quantum Lifshitz model. Finally, 
we conclude in section 6 with a conjecture that this suppression is universal for all Lifshitz 
theories in the absence of further constraining symmetries. 

2 The Green’s function in a scale invariant theory 

In a translationally invariant theory, the retarded Green’s function is naturally written in 
momentum space as GR{u},k). Furthermore, unitarity and causality demand that Gr is 
analytic in the upper half of the complex w-plane. For a scale-invariant theory, the conditions 
on the Green’s function are much stronger. In particular, for Lifshitz scaling symmetry with 
dynamical exponent z 

X —>• Ax, t —)■ A^t, (2.1) 

scale and rotational invariance demand that Gr cannot depend on co and k separately, but 
must have the form 

fe) = where Lb=-^. (2.2) 

\k\^ 

Here v is the energy scaling dimension and G{u}) is analytic in the upper half tD plane. 

Non-relativistic scale invariance by itself does not further constrain the form of Q{uj). 
However, additional symmetries can fix it completely. For example, relativistic conformal 
invariance (for the case z = 1) constrains Gr ~ where = —k^k^ = uY — \k\^. This 

is equivalent to taking the function 

Gc¥T = A{1 - , (2.3) 

where H is a constant. Similarly, full Schrodinger symmetry (for z = 2) [2-4] requires 

^Sch = 4l(l - 2mw)2^ (2.4) 
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where A is again a constant, and m is the eigenvalue of the mass-operator of the Schrodinger 
algebra. 

While the relativistic and Schrodinger cases are the most extensively studied, we are 
mainly interested in exploring the features of the function ^(cD) for Lifshitz models without 
additional symmetries using holographic methods. In general, Q will depend on the details of 
the model. However, some universal properties can be deduced in both the small and large a) 
limits. For cD —)> 0, the only dimensionful quantity that remains is |A:|. Hence Gji must behave 
as or equivalently 

Q(cu —)• 0) ~ const. (2-5) 

On the other hand, when uj —>■ oo, the dependence on |A:| drops out, and we must have 

Q(uj —>■ oo) ~ (2.6) 

As can be seen from (2.3) and (2.4), the z = 1 and Schrodinger z = 2 cases both satisfy these 
properties. 

The retarded Green’s function is in general complex, and this ought to be kept in mind 
when considering the limiting behaviors given above. Of particular interest is the general 
behavior of the spectral function x(w, k) = 2ImGij(w, k). For large oo, the spectral function 
scales as X ~ oo'^’', consistent with (2.6), as well as the relativistic and Schrodinger cases. 
The small oo limit, on the other hand, is more subtle. While scaling symmetry demands 
X ~ 2|/cp^^ Im^(ci)), with lmQ{Ld) approaching a constant as w —)• 0, this constant is in fact 
zero for the z = 1 and z = 2 Schrodinger cases. Moreover, for these cases x is identically 
vanishing for a range of w near zero. However, this no longer needs to be the case in theories 
with Lifshitz scaling, but without additional symmetries. Nevertheless, as we have shown 
in [5], in the latter case the spectral function is at most exponentially small in the limit 
(i) —>• 0, at least in the two-derivative holographic theory. What we will show below is that 
this exponential suppression of x remains robust, even when higher derivative corrections are 
included, as long as the perturbative expansion is kept under control. 


3 Holographic Lifshitz models 


The Lifshitz symmetry (2.1) can be realized in a gravitational background given by the metric 

(3.1) 


j 2 —dt^ -|- dp"^ ^ dx"^ 
^^ d+2 = “2 + 
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The boundary of the bulk spacetime is located at p = 0, while the horizon is at p = oo. For 
simplicity, we examine the scalar Green’s function, which can be holographically computed 
from the action of a bulk scalar (j){t,x,p). 

At the two-derivative level, the minimally coupled equation of motion for (p is simply 
(□ — m?)({) = 0. This system has been extensively studied, and the holographic computation 
of the retarded Green’s function is by now standard [14]. Working in momentum space and 
taking 

(/.(t, f, p) = (3.2) 


we find that ip{p) satisfies the Schrodinger-like equation —ip" + Uo'ip = 0 where 




p- 


P 


p.-2lz 


(3.3) 


and 


V = \ l + 


d + z 
2z 


(3.4) 


We can highlight the scaling properties of the solution by defining the dimensionless coordinate 


P = pW- 


(3.5) 


The Schrodinger-like equation now takes the form 

-iIj”{p)+ Uo{p)ip{p) = 0, ^o(p) = (3.6) 

In order to apply the AdS/CFT prescription for calculating the retarded Green’s function, 
we need to examine the solution near the boundary at p = 0 and as it approaches the horizon 
at p = oo. In the limit p —)• 0, the Schrodinger potential is dominated by the — l/4)/p^ 
term, and we find the boundary behavior 

z/;(p ^ 0) ~ Ap5-^ + Sp5+^. (3.7) 


Here we have used the convention that B is the coefficient of the normalizable mode, while A 
is the coefficient of the non-normalizable mode. For z > 1, Uq approaches —at the horizon, 
so the solution is oscillatory: 

'ip{p —)> oo) ~ 
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(3.8) 








For the retarded Green’s function, we take infalling boundary conditions, which correspond 
to setting 6 = 0. In this case, we find 


Giu) 


B 

A 


b=0 


(3.9) 


where the relation between {A, B} at the boundary and {a, 6} at the horizon is obtained by 
solving the Schrodinger problem (3.6). 


3.1 Bulk higher derivatives 

At the two-derivative level, the solution for ^(w) has been extensively studied, and analytic 
results may be obtained for z = 1 and z = 2 [1, 14]. However, as we emphasized in section 
2, scaling symmetry by itself does not fully constrain the form of the Green’s function. This 
raises the question of where the freedom of arbitrarily choosing the function G arises in the 
holographic dual. If we work within general relativity, there are two natural possibilities; the 
first is the choice of background metric, and the second is the form of the scalar equation. 
However, the metric (3.1) is essentially unique (up to coordinate transformations) once we 
have imposed Lifshitz scaling. This leaves us with modification of the equation of motion. 

From a bulk effective field theory point of view, it is possible to include higher derivative 
terms in the scalar equation. In momentum space, non-radial derivatives in the effective 
action show up as powers of u and k, while additional p derivatives lead to a higher order 
differential equation for ip{p). If there are no additional p derivatives, then the momentum 
space equation remains second order and can be brought into Schrodinger form just as above. 
This time, however, the effective Schrodinger potential in (3.6) generalizes to 

U{p) = ^ .2^^ + -' 2 ^/^ 

pi pZ-Z/Z pZ 

where the function / encodes the presence of the higher derivative terms. 

In principle, the procedure for extracting the holographic Green’s function is unchanged 
from the prescription of (3.9). However, the higher derivative terms affect the shape of the 
potential, and hence may change the boundary and horizon asymptotics and possibly also in¬ 
troduce additional classical turning points in the bulk. In order to get a better understanding 
of the asymptotics, we write out the expansion 

4/(^/5,p'/^)= E (3.11) 

^ id 

i+j>2 
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where i and j count the number of temporal and spatial derivatives, respectively. The re¬ 
striction i + j > 2 ensures that / only comprises the higher derivative contributions. Note 
that the coefficients Xij are dimensionless, although (after restoring units) we typically ex¬ 
pect Ajj ~ {ij, where i is some microscopic scale and L is the curvature scale of the 
Lifshitz bulk, such that £ <C T. 

Focusing first on the boundary at p = 0, we see that the behavior of the potential (3.10) 
remains dominated by the 1/ term, since i + j > 2 in the derivative expansion. Thus 
the boundary scaling behavior remains unchanged from (3.7), and the relation of the scaling 
dimension to zz is unaffected by the higher order terms. 

The horizon behavior, on the other hand, is considerably different. Since the horizon 
is located at p —)■ oo, and the expansion (3.11) in general contains positive powers of p, 
the successive higher derivative terms will become more and more dominant at the horizon. 
Furthermore, the potential will generically go to Too at the horizon, depending on the sign 
of Ajj of the dominant term. As a result, strictly speaking, the perturbative expansion of 
the scalar equation breaks down near the horizon. Nevertheless, we now argue that the 
holographic Green’s function can be extracted from the solution of the higher derivative 
equation in a controlled manner. 

3.2 Consistency of the higher derivative expansion 

At the two-derivative level, the Schrodinger potential (3.6) is monotonically decreasing as we 
move into the interior of the bulk, and there is a single classical turning point located at po 
where ?7o(po) = 0. For p < po, the solution connects to the power-law behavior (3.7) at the 
boundary, while for p > po, the solution is oscillatory, and infalling boundary conditions are 
chosen at the horizon. 

Ignoring the shift of in (3.6), there are two competing power laws in U, namely j 
and 1!, and the behavior of the solution depends on which of the power laws dominates 
at the classical turning point. We define the crossover point as p* = which is the location 
where the two terms become comparable. There are two distinct cases to consider: 

1. For Cj 3> the classical turning point is located at po ~ vjCj <C p*. This point 

is close to the boundary, and the 1/p^ potential ensures a power law behavior without 
exponential suppression. The holographic Green’s function is “featureless”, and behaves 
as ^ ~ . 
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2. For Cj <C the classical turning point is instead located at po ~ ^ p*. 

The Green’s function now probes deep into the bulk, and can have non-trivial features. 
Note that the wavefunction has exponential behavior in the region p* < p < poj leading 
to an effective decoupling of the boundary from the horizon [5]. 

We now consider the effect of the higher derivative terms, encoded in the function / in 
(3.11). Although this function dominates at the horizon, we nevertheless consider a formal 
perturbative expansion of the Schrodinger problem in the couplings Ajj. Of course, the higher 
order terms will dominate the wavefunction near the horizon. However, it is important to 
realize that the holographic Green’s function is not determined by the wavefunction at the 
horizon, but by its asymptotic behavior at the boundary. Infalling boundary conditions are 
needed at the horizon, but this can be imposed consistently at each order in the perturbative 
expansion. These infalling conditions will be seen in the boundary Green’s function, but will 
not dominate over lower orders in the expansion. 

Although a formal perturbative expansion can be used to solve the bulk scalar equation, 
the expansion of Q in the couplings \ij will only be sensible if the corrections can be kept small. 
Obviously this cannot be true globally, as the higher derivative terms typically dominate near 
the horizon. However, as one can see for example by using the WKB approximation (see 
section 4 and [5]), the holographic Green’s function only gives us information about physics 
between the boundary and the classical turning point po) where the wavefunction changes 
from exponential to oscillating behavior. Hence all that is necessary is to ensure that / 
remains small compared to the leading order potential Uq only for p < po- The specifics of 
this condition depend on whether we are in the high or low frequency regime. We consider 
these two cases separately: 

1. In the high frequency regime (w ;§> the dominant term in Uq is jf?. Since this 

term is decaying, while at the same time / becomes more important as we move away 
from the boundary, we only need to demand that / is small compared to /fP' at the 
classical turning point. This gives rise to the condition f{uipQ,pQ^^) <C P, which may 
be satisfied by taking <C 1, where we have assumed the expansion (3.11) along 

with the behavior of the couplings Ajj ~ (I'/L)*"’"-^"^. As we may see from (3.4), the 
scale of V is set by mL. Therefore, the condition for a valid expansion is equivalent to 
demanding m(. <C 1. We conclude that in this case, higher derivative corrections are 
under perturbative control provided the bulk couplings satisfy ml <C 1. This behavior 
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is very much like the relativistic z = 1 case, since in both situations the jfP' potential 

dominates up to the classical turning point. 


2. In the low frequency regime (tD <C we need to compare / with the term 

in Uq. Once again, we only need to consider the magnitude of / at the classical turning 
point. The condition is now /(capo, ^ ■> which gives rise to the requirement 

. (3.12) 


As oj is taken smaller and smaller, we need to take higher and higher order corrections 
into account. As a result, the perturbative expansion breaks down at small w, and 
results computed in this regime will not be robust against higher derivative corrections. 
Physically, what happens is that as ca —>■ 0, we probe closer and closer to the horizon, 
and it is precisely there where the higher derivative corrections dominate. 


Hence, as long as the scale of the bulk higher derivative corrections satisfies m^ <C 1, the 
perturbative expansion of the boundary Green’s function makes sense for dimensionless fre¬ 
quencies ca S> ■ For lower frequencies, the higher derivative terms start dominating. 

This feature of higher derivative terms becoming more pronounced at the horizon is not 
restricted to the Lifshitz background, but is in fact fairly general and shows up in, e.p., 
the pure AdS and Schwarzschild-AdS cases. While the pure AdS case tends to be robust 
against higher derivatives because of conformal invariance, more care may be needed in the 
case of holography at non-zero temperature [14-21]. Transport coefficients, such as the shear 
viscosity, may be extracted using the Kubo formula, which is evaluated at \k\ = 0 before 
sending a; —>• 0. Since this is consistent with (3.12), the perturbative expansion for transport 
coefficients is valid. At the same time, however, more care may be needed when analyzing 
general hydrodynamic modes, which are defined for both a; and k small, but nonzero (see e.g. 
[ 22 ]). 


4 WKB analysis of the spectral function 

In this section, we study the holographic spectral function of a probe scalar in Lifshitz space- 
time, in the presence of higher derivative corrections. To determine the effect of higher 
derivatives on the retarded Green’s function, we consider a probe scalar with an effective 
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potential of the form 


U = 


-1/4 


+ 


52-2/2 


-LO^ + 


E- 

i+j>2 




(4.1) 


The last term encodes an infinite set of higher derivative corrections to the equation of motion, 
where the (i,j) term corresponds to i temporal and j spatial derivatives. The size of the 
coefficients is expected to be set by a microscopic length scale I, so that (after restoring units of 
L) Xij ~ {ij. Since it is in general not possible to solve the corresponding Schrodinger 
equation for the potential (4.1) analytically, we will make use of the WKB approximation to 
obtain an approximate solution. This method can be used to calculate the imaginary part of 
the retarded Green’s function, which is proportional to the spectral function. After switching 
to the p coordinates defined in (3.5), the spectral function can be approximated by^ [5] 


K-hxaGR{uj,k) « |fe|2^Mime-2^e-2®. (4.2) 

e ->0 

Here A is a normalization constant and 

■s = [" 'iP]l0(p) + ^- (4.3) 

The additional l/p^ term is equivalent to an effective shift —)• + |, which is necessary 

for consistency of the WKB approximation for 1 potentials [7] . The integral is taken from 
a UV cutoff e to the classical turning point pQ. The WKB approximation for the imaginary 
part of the rescaled Green’s function defined in (2.2) is given by 


K-Hmg(cd) « lim (4.4) 

e ->0 

This expression is valid for a potential with only one classical turning point, such that the 
wavefunction is oscillating near the horizon and tunnels towards the boundary. Close to the 
boundary, the part of the potential leads to a power-law scaling of the wavefunction, 
which is stripped off by the factor of e~^'^ in (4.4). We can use (4.4) to determine the imprint 
of higher derivative corrections on the spectral function, provided that \ij < 0. In this case, 
the potential goes to —00 at the horizon, but the wavefunction still remains oscillating and 
we can consistently impose infalling boundary conditions. Later we will argue that (4.4) can 
in fact be used to provide a formal expansion for corrections with arbitrary sign. 

'^The additional prefactor of arises from letting e —^ which is the proper UV cutoff needed to 

cancel the log-divergence of the integral. 
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In order to perform a perturbative expansion of the WKB integral (4.3) in terms of Ajj, 
we need the higher derivative corrections to be subdominant compared to the other terms in 
U, at least in the domain of integration. We therefore demand 

(4.5) 


for all 0 < /> < po (see also the discussion in section 3.2). We can already see that this imposes 
an (h-dependent condition on the coefficients Xij, which we will make more explicit in what 
follows. 

We can now determine the leading order correction to Im Q (cu) by formally expanding the 
WKB integral in terms of the Ajj. At leading order, the higher derivative contributions are 
linear, so for our purposes it will be enough to drop the sum in (4.1) and only consider the 
effect of a single correction term with fixed (i, j). In a realistic model with a tower of higher 
derivative corrections, one may obtain a perturbative expansion for \mQ{Cj) by summing up 
the individual contributions, keeping in mind that if there are corrections at different order 
(e.g. a' and one may have to go beyond linear order to study the effect of all correction 

terms. 

A consistent expansion in Ajj- requires expanding both the integrand and the upper 
bound po, since the location of the turning point depends on the details of the correction 
terms. Writing S = + 5S, where is the two-derivative integral with Ajj = 0, we find 

(see appendix A for a rigorous derivation): 




p2 p2-2lz 




5S 




dp- 




(4.6) 

(4.7) 


where is the turning point for the case Ajj = 0, i.e. the solution of 


V 1 .2 


a2 “*"A2-2/^ 
Po Pq 


- = 0 , 


(4.8) 


and we expanded up to linear order in Xij. The large and small w-behavior of the unperturbed 


12 ~ 








integral was computed in [7]: 






5^°) (w < u^~^) ^ -zv + zv log ( 2 j/) + v{z - 1 ) log z + 




2 ( 2 - 1 );_ 1 

W 2-1 . 


zV 


( 2 ( 2 - 1 )) 


(4.9) 

(4.10) 


Let us now calculate the leading correction (4.7) in the same limits. For a) )?> , the 

unperturbed turning point lies at ~ which is well within the region where the l/p^ 
term dominates over Hence we can approximate the integral as 

(55; 


dp- 


2j^-ch2 


(4.11) 


Letting x = ujp/v, we find 


5S « 


1 - 2 \ t 


UJ 


/ dx- 

I ^ r\ 

u ^ 




(4.12) 


- 1 


For Cj'^ correction terms with j 7 ^ 0 are highly suppressed. After taking the UV cutoff 
e to zero, we therefore have 


<55: 




/I/2-1 


where 


^i—1 


a = / dx- 


^/^TT 


(4.13) 


(4.14) 


JO '2^/T^ 4r(^)' 

Using (4.13) and the unperturbed result (4.9), we arrive at the final answer 

A“^Im^(w ~ C = ( 2 z/)“^^exp \2u (l — + ''')] ’ (4-15) 

where the ellipsis indicates terms that are higher order in A. The scaling of Q with 
reflects the fact that at large frequencies, the Green’s function Gji{uj, k) = \k\‘^^^Q{Cj) becomes 
independent of k (see the discussion in section 2). The higher derivatives simply renormalize 
the numerical prefactor in a controlled way. The size of the higher derivative corrections at 
large Cj is controlled by 




i-2 


{m£) 


i-2 


(4.16) 
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where z > 2 is the number of temporal derivatives. Note that <C 1 is precisely what 

is required for the higher derivative corrections to be small up to the classical turning point 
PO) in the limit of large uj, as one can see by evaluating (4.5) at po in this limit, and noting 
that the unperturbed potential is monotonically decreasing. 

We now turn to calculating the higher derivative corrections in the case of small frequen¬ 
cies (w <C In this case, the unperturbed classical turning point lies at ~ 

We can split up the integral (4.7) in the following way [7]: Let p* = be the crossover 
scale, defined in the beginning of section 3.2, at which the two different terms in the potential 
(3.10), !j? and become comparable. Since a) <C we can then introduce a 

regulator scale pr such that p^ pr Pq^\ and split up the WKB integral in (4.7) as 

rpr 


<^Po 


rPo 


+ 


(4.17) 


' pr 


The first of the integrals above is taken over e < p < pr ^ Pq^\ so we can approximate the 
potential in this region as 

'2 1 
p- p^ 

On the other hand, the second integral is taken over Pr < p < Pq^\ so in this region we can 
write 


fj ~ -_^_-_L \ 

^ 52 + A2-2/Z ^ P 


(4.18) 


U 


1 


P 

Using these approximations, we find 
5S = 5Si + 5S2, 


,2-2/2 


-t52 + Ai>y+^'/^-2. 


(4.19) 




v 2 , 1 

7^“r^2V2 


dp- 

pr 2 


' p 2 — 2 /z 

Letting u = ^ the first integral can be written as 


1 -t52 


(4.20) 




4 
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vr+ 
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(4.21) 


where the integration bounds are —)■ 0 and Ur = {pr/p*)^^^ 1- In the small 

frequency limit uj <C , the correction term is highly suppressed unless z = 0. Hence we 
have 


<55i <5i,oyAoj^2i-2 
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J Uf 


i-l 

du , = + O 


VT+ 
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UJ 

, 1-2 


(4.22) 
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The remaining integral is divergent as —>■ oo. However, one can show that the contribution 
of the upper bound cancels with that from the lower bound of 882 , since Ur is after all a 
fictitious regulator scale. Hence the only contribution of (4.22) to 5S is due to evaluating the 
integral at the lower bound —)■ 0: 


^ + O 



(4.23) 


where 

/■“=o 7x5-1 

dj = / du =. (4.24) 

This contribution is hnite; in particular there are no loge terms, which would affect the 
boundary scaling. Similar to the large cu case, higher derivative corrections are controlled 
by terms of order ~ where n counts the number of derivatives. This becomes quali¬ 

tatively different when considering 6 S 2 , which captures the contribution of higher derivative 
corrections deep in the bulk. Letting x = we obtain 


6S2 = 1:L0 —^eijXijLO 2 )^ 


(4.25) 


where 


r 

J Or / 


J-1+ 


j-i 


ea = I dx 


(4.26) 


„ 2 -- 


' Pr/P q \ / \ _ ^ 

When expanding the lower bound in powers of pr/p^\ each term is designed to cancel with 
the corresponding contribution from 5Si. Instead of carrying out this cancellation explicitly, 
we can therefore let 

VttT 






dx- 


2 (^- 1 ) 


(4.27) 


'» P-;)r( "7(L7T )’ 

together with the prescription (4.23). Using (4.23), (4.25) and the zeroth order result (4.10), 
we arrive at the final answer 


K <C y 


l-z\ 


D exp 


—Co ^-^^E^Co) 


(4.28) 


where 




exp 


E{Co) = 




2 zy{l + difidjXojy^ ^ H-) 


zT 


2 ( 2 - 1 ) 


+ CjjAijO) ^ 2) _j_ . . . ^ 


(4.29) 
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Here the ellipses indicate terms that are higher order in A. We see that the higher derivative 
terms have two distinct effects: First, corrections with i = 0, which correspond to purely 
spatial derivatives, affect the overall normalization D of the spectral function. Second, and 
more importantly, higher derivative corrections with any i and j change the behavior of the 
spectral function as w —>■ 0, encoded in E{Cj). The w-dependent correction terms become more 
and more important at small frequencies, and eventually the perturbative expansion breaks 
down. This was to be expected, since at small the spectral function probes deep into the 
bulk, where higher derivatives dominate. However, recall that the coupling constants Ajj- are 
generically given by a ratio of a microscopic versus macroscopic length scale. A* j ~ (£/L)*+J-2. 
It is thus possible to keep the corrections in (4.28) small by demanding 

> w > . (4.30) 

This is precisely the bound we argued for in section 3.2. Since the condition (4.5), eval¬ 
uated at large a), also guarantees that iv/L <C 1 (see the discussion around (4.16) and 
section 3.2), there is a wide range of frequencies that satisfy the inequality (4.30). For fre¬ 
quencies within this range, (4.28) is a universal result: The spectral function behaves as 
~ exp (—const. • , and there are both constant and w-dependent corrections that 

can be computed order by order in perturbation theory. The naive limit tD —)• 0 is non- 
universal, since higher derivative corrections cannot be kept under control. 

The procedure for calculating higher derivative corrections to the spectral function out¬ 
lined in this section can in principle be applied to arbitrary corrections of the form (4.1). 
Note, however, that since we generally expect an infinite number of such corrections, going 
beyond leading order in ^/L may require expanding (4.3) to the appropriate order. 

Finally, let us comment on the sign of Ajj-. In the analysis above, we assumed that 
Ajj < 0, so that the wavefunction is always oscillating at the horizon, and no additional 
turning points are introduced. If Ajj is positive, the wavefunction is tunneling in the deep 
IR, leading to another tunneling contribution S'ir to the spectral function (4.2). At large 
enough p, the higher derivative corrections will always dominate the potential, so does 
not have a perturbative expansion in Ajj. This is simply a consequence of the fact that the 
potential in the IR is always sensitive to all of the (in principle infinitely many) coefficients 
that appear in the series of higher derivative corrections, and we cannot solve the equation of 
motion perturbatively in the IR. It therefore seems that one cannot trust our analysis in the 
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case of generic corrections with arbitrary sign. However, one can circumvent this problem in 
the following way: For any given tD, we can define a regulator surface at p = Phi<^) ^ Po, 
such that higher derivative corrections are still small at ph, i.e. 




Ph ^ Ph ■ 


(4.31) 


This guarantees that the wavefunction is still oscillating at ph, even though eventually higher 
derivatives may cause the potential to bend upwards again. Ignoring the (unknown) behavior 
of the wavefunction in the deep IR, we only impose infalling boundary conditions at ph, instead 
of /5 —7- oo. The surface at ph thus becomes an “effective horizon”, where the wavefunction is 
infalling: 

V^(p Ph) ~ (4.32) 

Here <h is an increasing function of p. The retarded Green’s function can then be computed 
using the usual formula (3.9), and the spectral function can be calculated approximately 
using the WKB-formula (4.2). On a practical level, this regularization prescription amounts 
to simply taking (4.2) for granted, and formally expanding the WKB integral S in Ajj, without 
worrying about the dynamics close to the horizon. 


5 Field theory models with z = 2 

As we have demonstrated, a holographic computation of the spectral function yields the 
universal low-frequency behavior x ~ exp(—const. • provided w is in the range 

(4.30), where the higher derivative corrections are controlled. From a field theory point of 
view, such an exponential behavior is not expected to arise at any finite perturbative order, 
but can show up non-perturbatively. This, of course, fits the framework of non-relativistic 
holography, where the field theory dual is expected to involve strong correlations. 

In this section, we explore two field theoretic models exhibiting z = 2 Lifshitz scaling. 
The first model is the quadratic band crossing model of [12], and the second is the quantum 
Lifshitz model [13]. Our strategy will be to identify phase-space regions with nonzero decay 
rates for bosonic quasi-particles, which, according to the optical theorem, will contribute to 
the imaginary part of the corresponding bosonic Green’s functions, and hence the spectral 
function. In both models, we confirm the presence of exponential suppression in the spectral 
function at small tD, in agreement with the holographic computation. 


- 17 - 


5.1 The quadratic band crossing model 


To set up the quadratic band crossing model, let us start with a massless Dirac theory in 
2 + 1 dimensions, with action^ 


S 


dxdt 


{i'yodo - ijidx - ij2dy) ^ 


(5.1) 


Here 'h = ('!/’i,V’ 2 )^ is a two-component spinor and 'I' = 'I'iyo- The 2 + 1 dimensional Dirac 
matrices are given by 



The interaction term in (5.1) is the only four-fermi term allowed for a two-component spinor. 
In the IR, this term is also the most relevant interaction term in the RG sense. At the 
Gaussian fixed point, this theory is conformally invariant with dynamical critical exponent 
z = 1. By setting the speed of light to unity, the theory contains only one control parameter, 
which is the interaction strength g. 

The quadratic band crossing model generalizes the above to form a scaling invariant 
model with z = 2 [12]. It does so by replacing the derivatives in the Dirac theory (5.1) by 
the following operators: 


ido ido + toV^, 

id^ 

idy -2t2dxdy, (5.3) 

where + dy and Iq, ti and t 2 are real parameters. After this substitution, we obtain 

a model with z = 2: 

S = j dxdt [70 [ido + toV^) + 7 R 1 {dl - dy) + 2 -f2t2dxdy] • (5.4) 

This action bears some resemblance with the original Dirac theory. However, in direct contrast 
to the Dirac theory, whose action only contains hrst-order derivatives, this model has a hrst 
order time derivative and second order spatial derivatives. As a result, space and time have 
different scaling dimensions, and it is straightforward to show that dimensionally [t] = 2 [x], 
^Note that we use signature (+, —, —) for the field theory. 
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corresponding to z = 2 at the Gaussian fixed point. In condensed matter systems, this 
model describes band touching points with quadratic dispersions, which have been observed in 
bilayer graphene (see for example the review articles [23-25]); a realization has been proposed 
in optical lattice systems, using high angular momentum orbitals [26]. 

Generically, the action (5.4) contains four control parameters: to, ti, t 2 , and the interac¬ 
tion strength g. However, we can set one of the three tiS to unity (say t 2 = 1) by rescaling. 
As shown in [12] and detailed in Appendix B, if we require SO(2) spatial rotational symmetry, 
then ti and t 2 must coincide. Furthermore, if a fermion particle-hole symmetry (i.e. charge 
conjugation) is enforced, then to must vanish. Here we will focus on the case with to = 0 and 
ti = t 2 = 1, which preserves both the spatial rotational and charge-conjugation symmetries. 
In this case, the action reduces to 

S = J dxdt [7oi<9o + 7i - dy) + 2 j 2 dxdy] 4' - , (5-5) 

and as shown in Appendix B, the free dispersion relation for the two bands is given simply 
by 

e±{k) = ±A:^. (5-6) 

It is worth emphasizing that most of our conclusions remain valid as long as |to| < |ti| and 
I to I < |t 2 |. As discussed in Appendix B, these inequalities ensure that the model has both 
particles and holes in the weak coupling limit (small g). 

5.1.1 Renormalization group analysis 

At tree level, the ■0^ term in the quadratic band crossing model is irrelevant (relevant) in the 
IR for systems above (below) 2-|-l dimensions. In 2-|-l dimensions, g is marginal at the tree 
level. A one-loop RG analysis indicates that a repulsive interaction (g > 0) is marginally 
relevant at IR, while an attractive interaction g < 0 is marginally irrelevant in the IR [12]. 
In the Dirac theory (5.1) on the other hand, the 0^ term is irrelevant (relevant) in the IR 
for systems above (below) H-l dimensions. In 1-1-1 dimensions, due to the special properties 
of the 1-|-1 conformal group, the 0^ term remains exactly marginal, before the system hits a 
Kosterlitz-Thouless transition. 

5.1.2 Boson correlation functions 

Although the model discussed above describes fermionic fields, bosonic modes can be con¬ 
structed from these fermionic degrees of freedom in the form of fermion bilinears. In the 
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(a) A one-loop contribution to the spectral func- (b) A live-loop contribution to the spectral func¬ 
tion, corresponding to the decay of a boson into tion, corresponding to the decay into three parti- 

one particle-hole pair. cles and three holes. 

Figure 1. Self-energy corrections for the boson modes. Here, solid lines represent fermionic propaga¬ 
tors and wiggly lines are boson propagators. 

particle-hole channel, we can build four different fermion bilinears (boson modes) 

hi = (5.7) 

with i = 0, 1, 2, and 3, and the fourth gamma matrix is given by 73 = ^ 707172 - Here 60 
is the fermion density operator and the other three bosonic operators can be used as order 
parameters for various symmetry breaking phases (nematic or quantum anomalous Hall) [12]. 
At the Gaussian fixed point, these bosonic modes have z = 2, which is inherited from the 
fermions. 

Additional bosonic modes can also be created in the particle-particle channel (e.g. ' 0 i' 02 )’ 
which are the order parameters for various superconducting states. In this section, we will 
only consider fermion bilinears in the particle-hole channel. These bosons can decay into 
particle-hole pairs and are thus expected to have a finite lifetime. Via the optical theorem, 
the existence of such decay channels is equivalent to a non-zero imaginary part of the two- 
point function (and thus the spectral function), generated by self-energy diagrams such as 
those shown in Fig. 1. 

Although it is challenging to analytically compute these diagrams, it is straightforward 
to prove that for a boson with momentum k, the imaginary part of each self-energy diagram 
can only arise when the energy a; of the boson is larger than a certain threshold. For each 
diagram, this threshold can be determined using energy-momentum conservation. 
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For example, the one loop diagram shown in Fig. 1(a) (the leading order correction) 
computes the scattering rate for a boson mode with energy ui and momentum k to decay into 
one particle with energy Up and momentum kp and one hole with energy w/j and momentum k^. 
Such a decay process can only take place when both the energy and momentum conservation 
laws are satisfied: 


k — kp kh^ 

2 ,2 

Lo — ujp oj}i — kp + kf^ ^ . 


(5.8) 


Here we used the quadratic dispersion relation (5.6). For fixed k, the momentum conservation 
law enforces a relation between the momentum of the particle kp and that of the hole kh, 
i.e. kp = k + kh- With this constraint, the energy of the particle-hole excitation kp + kj^ 
has a lower bound of k"^/2 (which is reached when kp = —kh = k/2). In other words, the 
energy conservation law can only be satisfied when w > A:^/2. As a result, for w > /c^/2, the 
boson can decay into a particle-hole pair, and thus have a finite lifetime, while for uj < A:^/ 2 , 
decay is kinematically forbidden. Thus, at the one-loop level, 0{g^), the imaginary part of 
the bosonic correlation function only arises for uj > /c^/2. This energy range is known as the 
particle-hole continuum. 

When the '^'^-interaction term is taken into consideration, the bosonic modes can decay 
through higher order processes (one example is shown in Fig. 1(b)). For these higher order 
diagrams, the same analysis can be utilized. At order the energy and momentum 

conservation laws imply that 


n+l 


n+1 


i=l i=l 

n+l n+l 

Y^ujp^-Y^ Uh, = + kl)> ^ . 


n+l 


k^ 


OJ = 


(5.9) 

*=i i=i *=i 

Here we consider the decay of a bosonic mode into n-|-l particles and n-|-l holes (see Fig. 1(b) 
for an example with n = 2). For fixed k, momentum conservation enforces a constraint on the 
momenta of the particles and holes. Given this constraint, the energy is minimized when the 
momenta are collinear, and the boson momentum k is equally distributed among the particles 
and holes. This results in a lower bound on the energy of kP'/2{n + 1). Thus the decay is 
kinematically forbidden unless u > k'^/2{n + 1). 
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This analysis demonstrates that, np to order of the imaginary part of the boson 

correlation fnnction only arises when the energy of the boson is above a threshold, uj > 
/c^/2(n + 1). Furthermore, this threshold goes down to zero for higher order diagrams as 
~ l/(n + 1). Thus, if we sum the diagrammatic expansion to infinite order (n —)• oo), we 
expect that the boson correlation function can pick up a nonzero imaginary part for any 
a; > 0. 

Finally, we are ready to extract the asymptotic form of the imaginary part of the self¬ 
energy correction at small oj. For iv k'^, the imaginary part can only arise via high order 
process where n ~ k'^j2uj. Therefore, we expect the imaginary part at energy lo and 

momentum k to scale as ~ g^^ ~ For sufficiently small g, this relation implies that 

the imaginary part of the self-energy correction decays to zero with the singular behavior 
~ g-const./o;^ where a) = iv/k'^ is the dimensionless energy. This matches the z = 2 low 
frequency behavior (4.28) obtained holographically. 


5.1.3 Dirac theory revisited and systems with higher 2 : 

We can repeat the kinematical analysis used above for similar models with arbitrary z > 1. 
In this case, for 0(g^^), the energy-momentum conservation law becomes 

n+1 n+1 

k = ^ kp, — ^ ., 

i=l i=l 

n+1 n+1 n+1 

^ + +) > 


k^ 


(2n -h 2) 


z-l ■ 


(5.10) 


i=l i=l 2=1 

For any z > 1, the lower bound for having a nonzero imaginary part depends on n, and goes 
to zero as n —)> 00 (i.e. when considering higher and higher order diagrams). Similar to the 
discussion above, after summing over all the diagrams to infinite order, we find that at small 
cj, the imaginary part of the self-energy scales as 


Imll ~ g 




(5.11) 


For small g, this indicates that Imll decays to zero as g-const.-o; where now to = 

cv/k^, in agreement with the holographic result (4.28). This suggests that the exponential 
suppression of the spectral function is a generic property of Lifshitz models at cj <C 

Note that for z = 1, the Dirac theory is recovered, and the fate of the system is funda¬ 
mentally different. As can be seen by substituting z = 1 into (5.10), the energy threshold 
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becomes independent of n. For any diagram, regardless of its order, the imaginary part arises 
only for u > k. After summing over all diagrams (to infinite order), the same lower bound 
of energy remains {uj > k). As a result, for z = 1 the imaginary part of the correlation 
function vanishes identically in a finite region ui < k, which is in sharp contrast to the z > 1 
case. This conclusion is consistent with a symmetry analysis, which tells us that at z = 1, the 
Lorentz and conformal symmetries require the bosonic correlation function to be proportional 
to (—a; + |A:|)“, where a is some scaling exponent. For non-integer a, {—uj + |fe|)“ is real for 
UJ < \k\, while the imaginary part arises for uj > \k\. For z > 1, however, the absence of the 
Lorentz and conformal symmetries allows for very different types of behavior. 

In summary, we find that models with z > 1 and z = 1 belong to fundamentally different 
universality classes. The case with z = 1 (i.e. Dirac) has been well understood with the 
help of conformal symmetry, which almost fully hxes the functional form of the correlation 
functions. However, for z > 1, the absence of conformal symmetry allows for richer structure 
in the correlation function. For arbitrary z > 1, we have presented an argument suggesting a 
characteristic exponential behavior g-const./aji/(^ £qj, imaginary part of the self-energy 
correction at low energy. 

5.1.4 Limitations of the analysis 

An exponential fall-off ~ g-const./w^/^^ ££^g spectral function all the way down to a; —)■ 0 
would correspond to an essential singularity of the two-point function at the origin. However, 
it is worth noting that there are two limitations of the analysis presented above. First, 
because we only considered the decay of bosonic modes into n -|- 1 particle-hole pairs without 
taking into account the renormalization of the vertex function (i.e. the renormalization of 
the coupling constant g), the above analysis is not expected to give quantitatively accurate 
results in the extremely low (or high) energy limit. This is because, as discussed above, in 
2+1 dimensions, the coupling constant g is marginally relevant or irrelevant (depending on 
the sign of g). For the IR or UV limit, the flow of g cannot be ignored. However, because g is 
only marginally relevant or irrelevant, the flow of g is expected to be slow (i.e. logarithmic). 
Hence there may exist a range for uj (i.e. uj is small, but not too small) in which the RG flow 
of g may be weak enough to be ignored, so that the analysis above can produce a reasonable 
estimate for the scaling behavior of Im H. 

Second, in the context of QFT, the perturbation series in terms of Feynman diagrams 
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is typically expected to be an asymptotic series. This means that our kinematical argument 
using loop diagrams only captures the behavior of the imaginary part of the self-energy 
correctly up to some finite order where N is large but finite. In particular, this 

implies that the scaling Im 11 ~ ~ ig only valid for n < N and thus for ui above 

some cutoff w*(A^). 

Both of these points suggest that while the exponential suppression of the spectral func¬ 
tion is a generic feature in a finite region where cD is small, the behavior in the strict limit 
ti) —7- 0 is model-dependent. This is consistent with the observation in the gravity theory, 
where the would-be singular behavior of the two-point function may receive significant cor¬ 
rections at very small lo from model-dependent higher derivative terms. 


5.2 The quantum Lifshitz model 


We now turn to the quantum Lifshitz model [13], which is a z = 2 generalization of the 
Klein-Gordon theory. We start with the action 


S 


dxdt 




(V2$)^ 


(5.12) 


Similar to the <h^-model, positive (negative) m corresponds to the disordered (ordered) phase 
respectively. The quantum Lifshitz model focusses on the quantum critical point between 
these two phases (at m = 0), at which the system is scaling invariant. At tree level, the 
dimensions of the various quantities are 


M = 2, [fe] = i, m = b] = 6-d. 


(5.13) 


where d is the number of spatial dimensions. The fact that [w] = 2[k] implies z = 2. 

In the non-interacting regime {g = 0), the bosonic field <I> correspond to free bosons with 
quadratic dispersion relation 


LO = 


(5.14) 


and the (free) two-point correlation function is 

-w)) = 

Similar to the previous model, we consider the decay of a bosonic mode with energy lo 
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(a) The two-loop sun-set diagram contribution to (b) A six-loop contribution to the spectral func- 

the spectral function, corresponding to the decay tion, corresponding to the decay into seven bosons, 

of one boson into three bosons. 

Figure 2. Self-energy corrections in the quantum Lifshitz model. Here, solid lines represent bosonic 
propagators. 

and momentum k into 2n -|- 1 bosonic modes (with energy uji and momentum ki where i = 
1,..., 2n + 1) in the 2n-th order diagram 0{g‘^^). Energy and momentum conservation imply 

2n+l 

k= kj, 
i=l 

2n+l 2n+l ,2 

(5.16) 

2=1 2=1 

Once again, we find that the decay can only take place for oj above a threshold, which 
approaches zero as n goes to infinity. 

Using the same analysis, we see that the imaginary part of the correlation function is 
nonzero for any finite u, and at small u the imaginary part is ~ with n ~ k"^/2oj. As 
a result, the imaginary part scales as ~ , and we again recover the non-analytic z = 2 

behavior ~ g-const./cj ^ 
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6 Discussion 


In our previous work, we found that the spectral function for a minimally coupled scalar in a 
Lifshitz background was nonzero, but exponentially small, in the low-frequency regime to 
[5] . The analysis presented here shows that this behavior is a robust holographic prediction for 
field theories with Lifshitz symmetry, in the absence of further constraining symmetries. For 
the classes of higher derivative theories we study holographically, we generically find that the 
spectral function is suppressed in the low frequency region as y ~ exp(—const. • 
so long as a) ^ where £ is the length scale at which higher derivatives become 

important. 

On the field theory side, the Lifshitz scaling symmetry is a priori not expected to lead 
to a universal 2-point function, and perturbative calculations do not reveal any similarities 
either between different field theories with Lifshitz symmetry, or with the holographic theory. 
However, we were able to show that in both of the held theory models considered here, a simple 
kinematical argument involving energy-momentum conservation and a resummation of loop 
diagrams reveals a similar exponential suppression as predicted by holography. Furthermore, 
this exponential suppression is expected for any held theory containing the following three key 
features: The existence of particles and holes, an interaction that allows for decay channels, 
and a dispersion relation with z > 1 scaling symmetry. Therefore, we expect our conclusion 
to be generic and applicable to a wide range of systems (with z > 1) regardless of microscopic 
details, in agreement with the holographic prediction. 

Although in both the holographic and the held theory calculation, the exponential sup¬ 
pression is a robust feature of the spectral function for small w, the strict limit w —)■ 0 is 
non-universal in both cases. In the holographic calculation, the model-dependence enters 
through higher derivative terms, which introduce corrections whose size can be quantihed 
precisely (see equation (4.30)). However, the precise regime of validity of the held theory 
calculations is less clear. In both of the models considered here, the how of the coupling 
constant g can no longer be neglected when taking the exact limit tD —?■ 0. Instead of just 
being a simple exponential, the exact (nonperturbative) spectral function will therefore have 
a more complicated dependence on tD. Naively, one may expect a dependence of the form 

liji Q ^ g(Lu)‘^ \ ( 6 . 1 ) 

In 2-|-l dimensions, the coupling g is marginal, and we expect g to depend only weakly on ui, 
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so that the spectral function still shows an approximately exponential behavior. It would be 
interesting to further study the renormalization group flow of g to make a precise statement 
about the range of u for which this is the case. Along the same lines, in order to put a 
precise lower bound on a), it would be important to account for the fact that the perturbative 
expansion is in fact only an asymptotic series (see the discussion in section 5.1.4). 

In our field theory calculation, we found that Im^ ~ go that exponential 

suppression in fact only arises for g( <C 1. It is important to note that this is not in con¬ 
tradiction to AdS/CFT being a weak-strong coupling duality. The strong coupling nature 
of the field theory does not necessarily mean that the parameter g has to be chosen large, 
but rather that strong correlations (for example seen as long-range interactions) may emerge 
dynamically. This feature is familiar from the standard case of relativistic AdS/CFT, where 
it is not g'YM itself that is taken large, but rather the’t Hooft coupling 5 ym-^ ;:g> 1. In order 
to better understand the relation between strong/weak coupling on the field theory/gravity 
side in non-relativistic AdS/CFT, it would be desirable to develop a more precise version of 
the holographic dictionary for this case. 

Although we have chosen not to consider higher derivatives in the radial direction p be¬ 
yond second order, this is in fact not a true limitation of the perturbative analysis. Assuming 
we are only interested in solutions to the higher derivative equation that are perturbatively 
connected to the lowest order (i.e. the two-derivative) equation, we may always eliminate 
higher derivatives by substituting in the lower order equations. Consider, for example, the 
addition of a fourth order term to the Schrodinger-like equation (3.6) 

-+ U{p)'4){p) = (6.2) 

We now rewrite this as p" = Up — and take two derivatives to obtain = (Up)" — 

Xp^^\ Substituting this in the right-hand side of (6.2) and working only to linear order in A 
then reduces the equation to second order 

- P"{p) + U{p)Pip) - X{lJ{p)P{p))" = 0{X^). (6.3) 

While this equation is no longer in manifest Schrodinger form, it can be so transformed if 
desired. Thus our analysis is in fact applicable to this more general case as well. 

As we discussed at the end of section 4, the perturbative expansion of the spectral func¬ 
tion in terms of higher derivative coefficients Ajj- strictly speaking only makes sense if these 
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coefficients are chosen such that no additional turning points are introduced deep in the bulk. 
However, we argued that our formal perturbation series can still be used even in the case 
of higher derivatives with “wrong” sign, i.e. for the case where the effective potential bends 
upwards at large p. It would be interesting to determine if in a realistic theory, there are 
constraints on the signs of the coefficients \ij, for example due to bulk causality or unitarity. 
It would also be interesting to study string theory embeddings of Lifshitz spacetimes, where 
the coefficients of higher derivative corrections can be determined exactly, and calculate the 
corrections to holographic correlation functions. 
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A Pertnrbative expansion of the WKB integral 


We would like to obtain an approximate expression for the WKB integral (4.3) for a potential 
of the form 

U{p) = Uo{p)+6Uip), (A.l) 

where 

^0 = ^ + -2-2/z ~ ^ 

and 6U represents a small correction to the potential. To be precise, we assume that 6U is 
subdominant compared to the other terms in the potential for all p between the boundary 
and the classical turning point. To guarantee this, it is sufficient to demand that 

5U{p)<^Cj^ for 0 < p < pq. (A.3) 


We can then expand the turning point as follows: 

Poit) = P^o^ (1 + t H-), 


where 


<R7(p; 


( 0 )^ 


<C 1, 




(A.4) 


(A.5) 
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and is the turning point of the unperturbed potential, i.e. Uo{p^q^) = 0. The relative size 
of 6U at the unperturbed turning point is what controls the higher derivative expansion. The 
WKB integral (4.3) can be written as S' = + 6S, where 


6S = S- S(°) 



dp 


Uo{p) + 5U{p) 


I 




dpJUoip). 


(A.6) 


One could attempt to simply expand the above expression formally in t and 5U, and it turns 
out that this does indeed give the correct result (4.7). However, this approach is problematic, 
since Uq goes to zero at Pq^\ and thus the formal expansion parameter dl7 jUo blows up at 
this location. The solution is to split up the integrals in (A.6) in a way that the integrand 
always has a well-defined expansion in terms of 5U. To do this, we shift the first integral by 
rescaling x = pp^Q^/po{t), so that the upper bounds of both integrals are identical. We can 
then combine both terms to obtain 


SS ^ I dx\l- {I + 2t) + 6U (x) + f dxJuo{x) ^/l + V{x) - 1 , (A.7) 


X^ X 


2-2lz 


where 


^ 5U{x) + ^^t-2t 

V{x) = --, 

Uo{x) 


(A.8) 


and we expanded to linear order in t. The first term in (A.7) is due to the shift of the lower 
bound of the first integral in (A.6). Assuming that lima,_>o x^(517(x) = 0, this term evaluates 
to i 2 t after we send e —)• 0. To compute the second integral, notice that although Uo{x) itself 
blows up at the upper bound, the ratio V(x) remains finite everywhere. Moreover, it is small 
by assumption, so we can expand (A.7) in terms of V{x): 


6U(x) -r - 2 - 2/2 

dS ^ ut + I dx- ^ 


i: 


t-2t 


2^ 


(A.9) 


The integral over the terms linear in t exactly cancels the term, and we arrive at the final 
result: 


SS: 


^^^\jUix) 


'2VW0 


(A.IO) 


This is the first order correction to the WKB integral in the presence of a perturbation 5U. 
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B The free quadratic band crossing theory 


In the absence of interactions {g = 0), the action for the quadratic band crossing model, 
shown in (5.4), describes two species of free fermions with dispersion relations 


e±(k) 


to ± 


,litj+tl) + {tl-t 

2 ) cos 4:9k 

V 2 



(B.l) 


where e+(/c) and €-{k) are the energies for the two species of fermions at momentum k 
respectively. The angle 9k is the azimuthal angle of k. 

In contrast with the 2+1 dimensional Dirac theory, where the SO(2) spatial rotational 
symmetry arises automatically (even if the speed of light is different along x and y, the 
rotational symmetry can be obtained by rescaling), the z = 2 model here in general only 
preserves a four-fold rotational symmetry (see the cosine term in the dispersion relation). 
Continuous rotational symmetry is only recovered at ti = t 2 [12]. 

For |to| < |ii| and |to| < |^ 2 |, h is easy to realize that e+ > 0 and e_ < 0. Therefore, 
fermions with energy e+ and e_ are particles and holes respectively (or say particles and 
anti-particles). These particles and holes in general do not preserve the symmetry of charge 
conjugation, because e+ except when to = 0- Thus charge conjugation only becomes 

a symmetry at to = 0. 

If |to| is larger than both |ti| and |t 2 |, then e+ and e_ will have the same sign, so that 
they are both particles (or both holes). If |ti| ^ |t 2 | and |to| takes a value between |ti| and 
|t 2 |, the system is anisotropic and along certain directions we have a particle branch and a 
hole branch, but along certain other directions, we have two hole branches (or two particle 
branches). There is one special case with |to| = \ti\ = \t 2 \- Here, one of the fermions has zero 
energy for any momentum, and this is known as a flat band. For a flat band, higher order 
spatial derivatives become important and will in general lift the energy degeneracy. Because 
a flat band is typically associated with an infinite density of states, it is usually unstable 
(towards a certain symmetry breaking ground state) when interactions are introduced. 

In the main text, we focus on the case with to = 0 and ti = t 2 = 1, but most conclusions 
remains valid as long as |to| < |H| and |to| < 1^21- For this special case, the dispersion relation 
reduces to e±{k) = as given in (5.6). 


- 30 




References 


[1] S. Kachru, X. Liu and M. Mulligan, Gravity Duals of Lifshitz-like Fixed Points, Phys.Rev. D78 
(2008) 106005, [0808.1725], 

[2] D. Son, Toward an AdS/cold atoms correspondence: A Geometric realization of the Schrodinger 
symmetry, Phys.Rev. D78 (2008) 046003, [0804.3972]. 

[3] K. Balasubramanian and J. McGreevy, Gravity duals for non-relativistic GFTs, Phys.Rev.Lett. 
101 (2008) 061601, [0804.4053]. 

[4] A. Adams, K. Balasubramanian and J. McGreevy, Plot Spacetimes for Gold Atoms, JHEP 0811 
(2008) 059, [0807.1111]. 

[5] G. Keeler, G. Knodel and J. T. Liu, Hidden horizons in non-relativistic AdS/GFT, JHEP 1408 
(2014) 024, [1404.4877]. 

[6] S. Leichenauer and V. Rosenhaus, AdS black holes, the bulk-boundary dictionary, and smearing 
functions, Phys.Rev. D88 (2013) 026003, [1304.6821]. 

[7] G. Keeler, G. Knodel and J. T. Liu, What do non-relativistic GETs tell us about Lifshitz 
spacetimes?, JHEP 1401 (2014) 062, [1308.5689]. 

[8] S.-J. Rey and V. Rosenhaus, Scanning Tunneling Macroscopy, Black Holes, and AdS/GFT Bulk 
Locality, 1403.3943. 

[9] J. Hartong and N. A. Obers, Horava-Lifshitz Gravity From Dynamical Newton-Cartan 
Geometry, 1504.07461. 

[10] D. M. Holman and B. Rollier, Warped Gonformal Field Theory as Lower Spin Gravity, 

1411.0672. 

[11] T. Griffin, P. Horava and G. M. Melby-Thompson, Lifshitz Gravity for Lifshitz Holography, 
Phys.Rev.Lett. 110 (2013) 081602, [1211.4872]. 

[12] K. Sun, H. Yao, E. Fradkin and S. Kivelson, Topological insulators and nematic phases from 
spontaneous symmetry breaking in 2d fermi systems with a quadratic band crossing, Phys. Rev. 
Lett. 103 (Jul, 2009) 046811, [0905.0907]. 

[13] E. Ardonne, P. Fendley and E. Fradkin, Topological order and conformal quantum critical 
points. Annals Phys. 310 (2004) 493-551, [cond-mat/0311466]. 

[14] D. T. Son and A. O. Starinets, Minkowski space correlators in AdS / CFT correspondence: 
Recipe and applications, JHEP 0209 (2002) 042, [hep-th/0205051]. 

[15] G. Policastro, D. T. Son and A. O. Starinets, The Shear viscosity of strongly coupled N=4 
supersymmetric Yang-Mills plasma, Phys.Rev.Lett. 87 (2001) 081601, [hep-th/0104066]. 


- 31 - 


[16] G. Policastro, D. T. Son and A. O. Starinets, From AdS / CFT correspondence to 
hydrodynamics, JFIEP 0209 (2002) 043, [hep-th/0205052]. 

[17] G. Policastro, D. T. Son and A. O. Starinets, From AdS / CFT correspondence to 
hydrodynamics. 2. Sound waves, JFIEP 0212 (2002) 054, [hep-th/0210220]. 

[18] G. P. Herzog, The Hydrodynamics of M theory, JHEP 0212 (2002) 026, [hep-th/0210126]. 

[19] G. P. Herzog, The Sound of M theory, Phys.Rev. DOS (2003) 024013, [hep-th/0302086]. 

[20] P. Kovtun, D. T. Son and A. O. Starinets, Viscosity in strongly interacting quantum field 
theories from black hole physics, Phys.Rev.Lett. 94 (2005) 111601, [hep-th/0405231]. 

[21] A. Buchel, J. T. Liu and A. O. Starinets, Coupling constant dependence of the shear viscosity in 
N=4 supersymmetric Yang-Mills theory, Nucl.Phys. B707 (2005) 56-68, [h.ep-tli/0406264]. 

[22] D. T. Son and A. O. Starinets, Viscosity, Black Holes, and Quantum Field Theory, 

Ann.Rev.Nucl.Part.Sci. 57 (2007) 95-118, [0704.0240]. 

[23] J. Nilsson, A. H. Gastro Neto, F. Guinea and N. M. R. Peres, Electronic properties of bilayer 
and multilayer graphene, Phys. Rev. B 78 (2008) 045405, [0712.3259]. 

[24] A. H. Gastro Neto, F. Guinea, N. M. R. Peres, K. S. Novoselov and A. K. Geim, The electronic 
properties of graphene. Rev. Mod. Phys. 81 (2009) 109, [0709.1163]. 

[25] V. Kotov, B. Uchoa, V. Pereira, F. Guinea and A. Castro Neto, Electron-electron interactions 
in graphene: Current status and perspectives. Rev. Mod. Phys. 84 (Jul, 2012) 1067-1125, 
[1012.3484]. 

[26] K. Sun, W. V. Liu, A. Hemmerich and S. D. Sarma, Topological semimetal in a fermionic 
optical lattice. Nature Physics 8 (2012) 67-70, [1011.4301]. 


- 32 - 


